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Abstract. 

A method, recently devised to obtain analytical approximations to certain classes of 
integrals, is used in combination with the WKB expansion to derive accurate analytical 
expressions for the spectrum of quantum potentials. The accuracy of our results is 
verified by comparing them both with the literature on the subject and with the 
numerical results obtained with a Fortran code. As an application of the method 
that we propose, we consider the meson spectroscopy with various phenomenological 
potentials. 
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1. Introduction 

Recently, in Refs. PJE] a method has been developed that allows one to obtain analytical 
approximations with arbitrary precision to a large class of integrals appearing in many 
diverse physical problems. The method is based on a carefully chosen expansion, which 
has been shown to converge uniformly and gives excellent results already at first order. 
In Refs. PI |2] such method has been applied to a few problems in Classical Mechanics 
and General Relativity, allowing one to obtain simple analytical formulas that work 
even in the non-perturbative regime. In this paper we consider a further application of 
the method, showing that it is possible to obtain accurate analytical expressions for the 
spectrum of a quantum potential by using it in conjunction with the WKB expansion. 

WKB is a semi-classical approximation, whose accuracy improves for large quantum 
numbers, based on an expansion in h of the Schrodinger equation: it has been widely 
used in the past to get approximate evaluations of the spectra of quantum potentials 
and it is generally employed at the leading order in h, although its extension at higher 
orders has allowed one to recover, in special cases, the exact energy eigenvalues j^j. 
On the other hand, for a generic potential function V(x) the calculation of the WKB 
correction to a given order can be done only numerically and, in general, no analytical 
expression can be found for the WKB integrals and, as a consequence, for the spectrum. 

The WKB method has already been employed also in the specific case of 
the phenomenological quark- antiquark potentials, using both non-relativistic [3] and 
relativistic [H] kinematics, showing a remarkably good accuracy even for the lowest 
quantum numbers. However, the solution of the WKB equation, although simpler than 
the original Schrodinger equation, has been again done numerically and only in the limit 
of large quantum numbers it has been possible to obtain the spectrum in a closed form. 

In this paper, by applying the method of Ref. [I] to the calculation of the WKB 
integrals we show that one can obtain an accurate analytical approximation to the full 
spectrum of the potential (see also Ref. The accuracy of the approximation depends 
both on the order to which the WKB expansion is considered and on the order to which 
our method is applied. We apply the method to a few phenomenological quark-antiquark 
potentials that have been used in the literature to calculate the spectra of heavy mesons. 

In Sect. E] we discuss our method and in Sect. El we first apply it to the case of 
the non-relativistic Cornell potential [7j and then to the case of a linear potential with 
relativistic kinematics. Finally, in Sect. El we briefly summarize our results. 

2. Formalism 

The standard WKB approximation is obtained by solving for the energy the following 
equation: 




(1) 
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where x± are the classical turning points of the potential V(x) and E = V(x±) is the 
energy. 

Let us concentrate on the integral 



J = J dx^JE-V(x). (2) 

In the spirit of the Linear Delta Expansion jH] we interpolate the potential V(x) as 

V s (x) = V (x) + 5(V(x) - V (x)), (3) 

where Vq(x) is a suitably chosen potential that depends on one or more arbitrary 
parameters Ai,A 2 ,... that we collectively call A in the following. Vs(x) reduces to the 
full potential for 5 = 1. We want to perform an expansion in 5 without moving the 
inversion points: for this reason we impose that x± be the inversion points also of the 
potential V (x). As a result, the energy E that the particle would possess if it were 
moving only in the potential Vq(x) is given by E = V (x±). 
We therefore introduce a new integral 

J® = l X+ dx JE - V (x) [1 + 5A(x)} 1/2 , (4) 

J X- 

where 

Of course, one has J^ =1 "> = J. 

We then treat the term proportional to 5 as a perturbation and perform an 
expansion in 5, 

oo 

Jl*) = £ 5 n J n . (6) 

n=0 

Such expansion will converge uniformly when |A(x)| < 1 in the region x_ < x < x + PQ. 
This condition selects a particular region in the parameter space A; however, maximal 
convergence is achieved when the Principle of Minimal Sensitivity (PMS) ■9j is used, 
i. e. when, given the iV-th order approximation to J , 

N 

J(N) = J2^' ( 7 ) 

n=0 

the condition 

dJ (N) /d\ = (8) 

is enforced. Eq. (JSJ) fixes the parameters A, which are then inserted back in Eq. (|7j) 
yielding the approximant of J . 

Notice that if the potential V (x) is chosen appropriately it will be possible to 
calculate analytically each term in the expansion. 



Analytical approximations to the spectra of quark- antiquark potentials 

tOr 




Figure 1. Effective potential for a = k = L = fi = 1. 



3. Applications: meson spectroscopy 

3.1. N on-relativistic kinematics: the Cornell potential 
Following |H we write 
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The effective potential in this case is given by 

k L 2 

V eS (r) = h ar + - — - 

r 2[ir z 



(9) 



(10) 



where the last term is the centrifugal barrier and r± are the classical inversion points. 

As one can see from Fig. ^ this potential is not symmetrical with respect to the 
minimum, which we call tq. We therefore write 




E H ar — 

r 2[ir l 



f r + 






+ / dr* 
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and evaluate the two integrals separately. 

Following the procedure outlined in the previous section, we can write 



fro 

J±= T dr 

J r± 
rro 



\ 



2fi 



E H ar 

r 



L 2 
2^ 



f ° dr^/2/u [E 0± - V 0± (r)]^l + 5A ± (r) 

Jr± 



(12) 
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Figure 2. % error of J(\)+ + J(i\- fora = K = L = ^ = lasa function of E 



where 



and 



A+(r) 



E - E ± + ^ -ar 



2 fir 



£5 



0± 



^o±(r) 



^o±(r) = - r o) 2 ' ^o± = V 0± (r ± ), 



(13) 



(14) 



A± being arbitrary parameters. 

Provided that |A±(r)| < 1 for all r e (r±,ro), one can expand the expression for 
J± in powers of A±(r). At first order one obtains 



i± 



A 



where 



n 



i± 



± 



i± 



<X_/6 1_/ c 



(15) 



(16) 



^o±j ; ^i± an d Wi-t being analytic functions of r± and r given in the Appendix. 

By minimizing Jix\± with respect to the corresponding parameter \± one finds the 
optimal values, 



Ah 



I2H 



i± 



H 



o± 



(17) 



and, as a consequence, the optimal first order approximation to the integrals in Eq. (fTT| 
turns out to be 



J{i)± — \l'2n\l'rto±Hi±. (11 

In Fig. El one can see the typical accuracy one gets using this approximation. 
One has now to impose the WKB condition 



J m+ + J { iy = vr [n+ - 



(19) 
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and solve for E = E n . In order to simplify this equation we expand the left hand side 
in 1/E n to leading order: by inspecting the resulting expression one can easily check 
that the WKB condition can be satisfied by expanding the energy as a series in 1/n, 
i. e. assuming the following ansatz 



£„ = 5> fc Tl+1/2 
n=0 V 



L(8 + 17tt) 



2/3-/c/6 



(20) 



32\/27r J 

and determining the coefficients c& so that the quantization condition is obeyed. These 
coefficients can be recursively calculated using any computer albegra program: stopping 
at order 10 in the 1/n expansion one finds 

2 2 / 3 a 2 /W 
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(216) 
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(2 Id) 
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3k 



(8 - 2l7r)V7r - 2 - 48(tt - 2)^ 



-4a 
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(21$) 



c 9 



18432L 3 /V/6(2(7r - 2)) 7 /i%23/i2 r 3/2 
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-3v^ 
+448K 2 /i 2 7r 2 r 2 ) 



-256aL Vr (3Vtt - 2(4 - 3vr) + 2tt 3/2 
97tt))L 4 - 16K/i7r(777r - 24)r L 2 



(2 In) 



cio = 49152v ^ L 3 2 {( 8 + 7r )t 64 + ^(^^ " 144 )] L " + 32 ^ r o ^ 64 " 37r ( 32 + 37r )] 

-16a[16 + 7r(137r - 66)]r 2 ] L 2 + 256kVtt 2 (8 + 257r)r 2 } . (21*) 

In Table Q we compare the analytic spectrum of the cc and bb systems — obtained 
from Eqs. (|2Uj) and (|3.1|) — to the spectra obtained from the exact solutions of the 
Schroedinger and WKB equations, respectively jl]. 

As one can see from the table, in the case of the cc mesons the error in the analytic 
expressions is at most 1% for the lowest state and it is rapidly decreasing with the 
excitation energy; in the case of the bb mesons one goes from a 2.8% error in the ground 
state up to a 0.06% error for the 2P state. Note that the accuracy of the analytic 
approximation can be arbitraly improved simply by extending the order in the 1/E n 
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Table 1. Masses (in MeV) of the lowest cc and bb states from the exact solutions 
of the Schroedinger and WKB equations [3J and from the analytical approximation of 
Eqs. d2Hl> and (EH). Quark masses and potential parameters are taken from Ref. |10j . 



States 


Exact 


WKB oxact 


WKB ana i y tic 


cc 








IS 


3067 


3062 


3097 


2S 


3693 


3691 


3722 


IP 


3497 


3497 


3463 


2P 


3991 


3991 


3980 


ID 


3806 


3806 


3791 


2D 


4242 


4242 


4235 


bb 








IS 


9448 


9439 


9705 


2S 


10007 


10003 


10131 


IP 


9901 


9900 


9865 


2P 


10261 


10261 


10255 


ID 


10148 


10148 


10098 



and 1/n expansions and that similar analytic formulae can also be obtained for any 
observable by expanding its WKB expression. 

It is also quite straightforward to get simple formulae in the asymptotic regimes. 
For instance, if we assume n>L>l, the energy is approximated by our formula as 

a 2/3 

E n « c n 2/3 « 2.2245^n 2/3 , (22) 

which is extremely close to the result of Eq. (10) of Ref. j3]: 

a 2 / 3 , 

E n P3 2.2309— - K n 2lz . (23) 

On the other hand, in the regime L > n > 1 we have to make explicit the dependence 
of ro upon L and then expand around L = oo; we obtain: 

2/3 

E n « 1.62176^L 2 / 3 , (24) 

which again can be compared with the formula (9) of Ref. jl]: 

o 2 / 3 , 

E n » 1.5^L 2 / 3 . (25) 
3.2. Relativistic kinematics: linear potential 

A generalization of the Schrodinger equation to include relativistic effects can be done 
by starting with the Bethe-Salpeter equation and making a few approximations, namely 
assuming an instantaneous local potential and neglecting spin and the coupling of 
"large" to "small" components in the relativistic wave function. The resulting equation 
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also known as spinless Salpeter equation — has exactly the form of a Schrodinger 



equation with the kinetic term replaced with J — V 2 + m\ + y — V 2 + m\. 

Although the relativistic Schrodinger equation has been used in the literature 
with many realistic phenomenological potentials, in the following we shall consider for 
simplicity the linear potential V(r) = fi 2 r. In the relativistic case the WKB condition 
turns out to be slightly modified and it can be written as [3] 



where 



with 



J 



J 



( I 3 N 

*>+2 + z 



r+ 



dr 
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(E — fi 2 r) 2 m\ + m\ 
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In the limit mi = m2 
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2VT 



m, Eq. 1)26)1 can be integrated exactly to give [Hj 
E 



m log 



2m 



+ 



/ E 2 
Am 2 



(26) 



(27) 



(28) 



(29) 



Before dealing with the general case m\ ^ 777.2, we first want to show that Eq. ()29|) can 
be used to obtain an analytical approximation to the spectrum. We make the ansatz: 

/: = > < r, | •• - -I + d k [n + - + - ) log ( n + l - + 7 ) ^ . (30) 



^ I / / 3 
2^SCfc[n + - + - 

fc=0 



4 



The unknown coefficients and cffc in this expansion are determined by inserting 
Eq. flUJ) into Eq. (JH and by Taylor expanding in l/(n + 1/2 + 3/4) < 1. 



E„ 



n.l 



m 2 + /i 2 (3 + 21 + An) rr + m 2 log (^L) 
/iv / 3T2T+4n v /7r 
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m log ( ^ 
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m 2 + fi 2 (3 + 2/ + 4n) vr 
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2 log(fe) 



+ • • • • (31) 

/i 3 (3 + 2Z + 4n) 3/ V/ 2 

This formula has been obtained truncating the sum into Eq. ()3Uj) to k — 2. In Fig. |3] 
we compare the numerical solutions of Eq. (|26|) with the analytical approximation of 
Eq. (|3T|). using m = /i = 1 and Z = (solid curve). As one can see the accuracy is very 
good. 

We now want to extend our previous results to cases where the quantization integral 
cannot be performed analytically, such as for m x 7^ m 2 . Following the procedure outlined 
in Sect. EJwe can rewrite Eq. (|2*7|) as 



J(S) 



dr 



\ 



(E - X 2 r) 



+ 5 



(E — /i 2 r) 2 
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A(E - fx 2 ry 



(32) 
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Figure 3. Error over the spectrum denned as £ = (E° xact — fi^pproxy gexact usm g 
the analytical formulas obtained and taking m = jj, = 1 and I = 0. 



where A is an arbitrary parameter and 5 is used for power counting. For 5=1 the 
expression above reduces to the one in Eq. ()26j) . £?q is fixed by asking that the zero-th 
order term in 5 vanishes at the inversion point and its value is 

E = A 2 r+. (33) 

If we expand Eq. (I32J) in powers of 5 at order iV and then set 5 = 1, we obtain an 
expression that depends upon A. We again take the approximant J^n) and impose 
dJ~(N)/dX = to minimize such dependence. To leading order in S one obtains 



A 4 = /i 4 



-8 (Am) 2 (E + 2m) + E (E - 2m) (E + 6m) + 8(Am) 2 £ log(2£m) 



E(E - 2my 

where Am = (mi — m 2 )/2 and m = (m x + m 2 )/2. 
Thus, to order S we find 
E — 2m 



(34) 



J(X) 



A^Efi 2 



{E (E - 2m) (E + 6m) - 8(Am) 2 [E + 2m — E log(2)] 



1/2 



+8(Am) 2 E [log(.E) + log(m)] } ' " . (35) 

In Fig. [?] we compare the numerical result for obtained using mi = 1, m 2 = 2 and 
)U = 1, with the one of Eq. (}35|) . as a function of the energy E. Eq. (}35j) provides a quite 
simple and precise approximation to the exact integral. 
By introducing the ansatz 



l/2-fc/2 



/ Z 3 
+ d fc \n + - + - 



-l/2-fc/2 



log f n + ^ + | 



(36) 



into Eq. (|35|) and Taylor expanding the WKB condition in l/(n + //2 + 3/4) 1 one 
obtains to leading order 

6m 2 + 3/i 2 7r + 2//i 2 vr + 4// 2 n7r - 2(Am) 2 log(m/i) - (Am) 2 (-2 + log(167r)) 



E 



n. i 



fiy/3 + 21 + 4ny/7r 
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Numerical 

Analytical (first order) 




Figure 4. Comparison between the exact numerical value and the analytical 
approximation at first order (Eq. jSSJl) of the quantization integral ', taking mi = 1, 
77i2 = 2 and ji = 1. 

(Am) 2 (4mV3 + 21 + An -^(3 + 21 + An) 0r) log(| + f + n) 
H 3 h • • • • (37) 

/i 2 (3 + 2/ + 4n) 5 7r 

The dashed curve of Fig. |3] displays the error over the spectrum obtained using 
eq. (}3T|) . showing that accurate predictions can be obtained with this simple formula. 

4. Conclusions 

In this paper we have applied to the integrals occurring in the application of the WKB a 
simple method, which allows one to express them analytically and to arbitrary precision 
in terms of elementary functions (a formal proof of the convergence of such method has 
been previously given in Ref. PQ). This procedure is advantageous even in the case (not 
so frequent) in which the integral can be performed analytically, since the exact result 
will in general contain special functions, which are usualy more difficult to handle. 

Successively, the approximate WKB quantization condition obtained in the first 
stage has been inverted, thus providing analytical formulas for the spectrum of a particle 
in a given potential. 

In this paper we have considered two physical applications, which are relevant in 
meson spectroscopy, showing that in both cases excellent accuracy is obtained. The 
precision of our results essentially depends upon the intrinsic precision of the leading 
order WKB approximation, since the WKB integrals can be approximated analytically 
to arbitrary precision using our method. As a matter of fact we have also applied our 
method to one-dimensional anharmonic potentials (not discussed here) — using a higher 
order WKB calculation — and obtaining highly accurate formulas for the spectrum. 
Given the simplicity of our method we believe that it can be applied to a wide class of 
problems. 
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Appendix A. 



In this Appendix we provide the explicit expressions of the turning points and of the 
minimum of the Cornell potential (fTUj) and also the functions TLq± 1 TCf ± , 7~L b l± an d 7~C 
introduced in Sect. 13.11 

For the turning points we follow Ref. jl]: 
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;i-2) 



The functions 7i entering the first order approximation to the WKB integral of Eq. ()15j) 
are defined as follows: 
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